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Abstract 

Kontsevich's work on Airy matrix integrals has led to explicit results for 
the intersection numbers of the moduli space of curves. In this article we 
show that a duality between k-point functions on N x N matrices and N- 
point functions of k x k matrices, plus the replica method, familiar in the 
theory of disordered systems, allows one to recover Kontsevich's results on the 
intersection numbers, and to generalize them to other models. This provides 
an alternative and simple way to compute intersection numbers with one 
marked point, and leads also to some new results. 
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1 Introduction 



After Witten's celebrated conjectures [1, 2] on the relation between inter- 
section numbers on moduli spaces of curves and the KdV hierarchy, and 
Kontsevich's proof [3], the literature on the subject, both from the point of 
view of mathematics or from its string theory relationship, has become con- 
siderable. We want to add here a new method, duality plus replica, which 
allows one to recover easily some results which are difficult to obtain by 
fancier methods, and provides some new results as well. 

In a previous article [4] we have used explicit integral representations for 
the correlation functions [5, 6, 7] for a Gaussian unitary ensemble (GUE) 
of random matrices M in the presence of an external matrix source. The 
probabililty distribution for N x N Hermitian matrices is 

P^(M) = -Le-^*-^^'-^*--^^ (1) 

From this representation we have obtained the correlation functions of the 
'vertices' 

V{ki, kn) = ^< trM'^'trM''^ ■ ■ •trM'=" > (2) 

As is well-known from Wick's theorem (when the source A vanishes) such 
correlation functions are just the numbers of pairwise gluing of the legs of the 
vertex operators. The dual cells of these vertices are polygons, whose edges 
are pairwise glued. This generates orientable surfaces, discretized Riemann 
surfaces whose genus is related to the power of l/N"^. 

Okounkov and Pandharipande [8, 9, 10] have shown that the intersection 
numbers, computed by Kontsevich [3], may be obtained by taking a simul- 
taneous large N and large ki limit. In our previous work we have used the 
exact integral representation valid for finite N of those vertex correlation 
functions, and obtained explicitely the scaling region for large ki and large N 
by a simple saddle-point. This led to a practical way to compute intersection 
numbers from a pure Gaussian model, much simpler than Kontsevich's Airy 
matrix model. 

In this article we want to show that duality and replica may be used 
also to recover easily earlier results, to establish some new ones and give 
support to Witten's conjecture. In order to make this article self-contained 
we will add appendices in which we rederive some of the steps leading to the 
representation that we are using : 

(i) explicit formulae, valid for arbitrary N and arbitrary source matrix A for 
the average 

t/(si,---,Sfc) = (Wi^---IVe^^^) (3) 
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which rely on the Itzykson-Zuber formula [11, 12] (appendix A). 

(ii) a duality representation for the average of characteristic determinants 

(det(Ai-M)---det(Afe-M)) 

in terms of another GUE integral, but in which the random matrices are 
k X k , whereas the initial problem involved N x N matrices (appendix B) 
[13]. (This duality seems to be a simple reflection of the open string/closed 
string duality [14]). With the help of these two kinds of results we may 
proceed to the replica approach based on the simple relation 

n^O n OA A — D 

Therefore after applying these two steps we end up looking for an n goes to 
zero limit on matrix integrals whose size vanishes with n, a very different 
problem from the familiar large N limit. Its Feynman graph representation 
connects directly with Riemann surfaces with marked points of maximum 
genera. 



2 The duality plus replica strategy 

We first consider the average of products of characteristic polynomials, de- 
fined as 

1 ^ 

Ffe(Ai, Afc) = ^ < n det(A« - M) >a,m 

1 k 

^ _L fdMll det{Xi ■ I - M)e-Tt>^M2+NtrMA 

where M is an iV x Hermitian random matrix, A a given Hermitian ma- 
trix, whose eigenvalues are (ai, ■ ■ ■ , fln) and the normalization constant 
of the probability measure (for A = Q). We have shown earlier [13] that this 
correlation function has also a dual expression. This duality interchanges A", 
the size of the random matrix, with /c, the number of points in F^, as well 
as the matrix source A with the diagonal matrix A = diag(Ai, A^). This 
duality reads [13] 

FkiXu A.) = ^ / dB n[det(a, - tB)]e-^'^^^-'^^' (6) 

J j=i 

where A = diag(Ai, A^) and B is a k x k Hermitian matrix. (The nor- 
malization is now on GUE ensembles of Hermitian k x k matrices Zk — 
J d5exp(— |trB^)). The derivation is reproduced in appendix B. 
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If we specialize this formula to a source A equal to the unit matrix, a 
trivial shift for the original N xN matrices M, which has the effect of making 
the support of Wigner's semi-circle law, for the density of eigenvalues of M, 
to lie between and 2 rather than (-1,-1-1), the formula (6) involves 

det{l-zB)^ = exp[NtT\n{l - iB)] 

= exp[-iNtrB + ^tiB^ + i^trB^ ^ ] (7) 

The linear term in B in (7), combined with the linear term of the exponent 
of (6), shifts A by one. The B"^ terms in (6) cancel. In a scale in which the 
initials Afc are close to one, or more precisely N'^^^(Xk — 1) is finite, the large 
N asymptotics of (6) is given by matrices B of order A^~^/^. Then the higher 
terms in (7) are negligible and we are left with terms linear and cubic in the 
exponent, namely 

Ffe(Ai, Afe) = e>''' J dBe'T''^'+^^''^^^-'\ (8) 

So finally for this edge behavior problem, in which the matrix B is of order 
N~^^^, and A — 1 of order iV~^/^ we can rescale B and A — 1 to get rid 
of N. The result is nearly identical to the matrix Airy integral, namely 
Kontsevich's model [3], which gives the intersection numbers of moduli of 
curves. The original Kontsevich partition function was defined as 

Z^^J dBe-"'^''"+-s''^' (9) 

where Z' = J dBe"^^^^ . The shift B ^ B + iA, eliminates the B^ term and 
one recovers (8) up to a trivial rescaling. 

Let us apply this to a one-point function. In view of the replica limit we 
specialize those formulae to 

< [det(A - M)]" >A,M=< [det(l - iB)f >a,b (10) 

where B is an n x n random Hermitian matrix. For this edge problem we 
have chosen for the source matrix A, the N-dimensional unit matrix, whereas 
A is a multiple of the nxn identity matrix : A = diag(A, A). Note that the 
average in the l.h.s. of (10) is meant for the N x N GUE ensemble, whereas 
in the r.h.s. the average is performed on n x n Hermitian matrices with the 
weight (8) or (9). 

The strategy that we will used is thus to take the integral (9) for A = 
A X 1 and expand it in powers of the cubic term. The formulae that will be 
established, being exact for finite n, are easily continued in n and allow one 
to take the n — > limit. The method relies on explicit exact representations 
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of Gaussian averages [7, 15] in the presence of an external matrix source (in 
appendix A the main steps of the derivation have been recalled). As a result 
we will obtain formulae for quantities such as 

lim - < (trB'^)^ > 

computed with the (normalized) Gaussian weight Z~^exp(— ^trB^) onnx n 
Hermitian matrices ; thereby this procedure gives the values of the intersec- 
tion numbers of the moduli of curves for a number of cases. 



3 Evolution operators and replica limit 

We will rely on explicit expressions for 

C/A(si,---,Sfe) = -(tre^^''...tre^^'^) (11) 

Th 

for a probability measure on n x n Hermitian matrices in the presence of an 
Hermitian matrix source A, whose eigenvalues are ai, • • • , a„. The average is 
thus defined with the normalized weight 

Then one has (see appendix A) for the one-point function 

Ua{s) = i(tre^^) = — / — e-/^ A^u-a^ 
^ ^ ' ns 7 2i7r Y u - a„ ^ ^ ' 

the contour, in the complex u-plane, circling around the n poles a^^. We have 
chosen the normalization ^7^(0) = 1. This formula simplifies for a vanishing 
external source (although a non-vanishing source simplifies the derivation) 

ns J 2m u 
This representation leads to a simple continuation in n, and to an expansion 
in powers of n. For instance it gives 

e^^ rdu 

lim U{s) = — i — e^"/^ log (1 + -). (15) 

This last contour integral reduces to the integral of the discontinuity of the 
logarithm, giving readily 

lim U{s) = (16) 

\2\) 



and thus 

Note that in terms of Feynman diagrams with double hues, the hmit n 
selects the diagrams with one single internal index all along the lines of the 
diagrams. Those diagrams correspond to a surface of maximum genus for a 
given number of vertices. 

The same strategy works for higher point-functions. The k-point function 
is given (for a vanishing source) by 

t/(si,---,Sfc) = -(tre^^^---tre«^B) (18) 
n 



= (_l)fc(^-i)/2eEt ilfl ^e^>^'^/'^ n(l + det 

^ U/j, Ui ~\~ Si Uj 

Again the continuation to non-integer n is straightforward and leads to 
lim U(si, • • • , Sfc) = (-l)'=('=-i)/2eE^ i (19) 

n— >0 



k 



Ui + Si — u 



The calculation of the contour integrals is more cumbersome, but all the 
integration can be done explicitly to the end and give a remarkably compact 
result (for more details see appendix C). If we denote 

a = Si + --- + Sk (20) 

they lead to 

hm U{s,, . . . , 5fc) = - n 2 sinh ^. (21) 

This function is of course the generating function for the n = limit of 
^(trBPi • ■ ■ trBP'') by expanding in the Sj's. Selecting the coefficients of equal 
powers for every Si, for instance of (si ■ ■ ■ Skf', we find 

lim i < (trBY«- >= ^""^;:''y-^>' ('0-'] (22) 

all other powers of lim„^o ^ < (trB^)'^ > vanishing unless k — 2 (mod4). This 
leads to the intersection number of the moduli of curves with one marked 
point. Indeed , following Kontsevich, these numbers are given by 

-logZ = ^t;(r;) (23) 
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with 

ti ^ (-2)-(4W6 '-Q(2m + (24) 
^ 

from which the above result (22) provides 

<^^-^>^=(2^ (^ = 0,1,2,...). (25) 

These numbers agree with the values of the intersection numbers computed 
earlier by Kontscvich, Witten and others [3, 1, 16, 17]. 

Clearly the method allows one to compute more than that. For instance 
one can derive as well 

lim i < (trB^)*-' >= ^(4p - 3)! ( " ^ ) (26) 

These vertices appear in the higher Airy fuctions and they are related to the 
intersection numbers of Witten's top Chern class. [1, 3] 



4 Application to intersection numbers of top 
Chern class (p-spin curves) 

The partition function for the p-th generalization of Kontsevich's Airy matrix- 
model is defined as 

Z^ — [ c^Sexpf^— tr(BP+i - Ap+^) - tr(B - A)Ap] (27) 
Zq J p+l 

normalized by the Gaussian part of the integrand (after the shift B ^ B + A 
which cancels the linear terms in B of the exponent) : 

Zo- /ci5expEtr^AjBAP-j-iB] (28) 

The 'free energy', i.e. the logarithm of the partition function, is the gener- 
ating function of the generahzed intersection numbers (H'^mj) for moduli of 
curves with 'spin' j [3, 18] 

^=E<nvr>nf^ (29) 

dm,j nT-d ^,3 "^d' 

where [3] 



6 



, , OA rn—1 1 

tm,j = (-p)T^ n i^P + J + 1)^^ A^^^ (30) 
According to Witten [2] the intersection numbers is given by 

< W-Tn.„J„>.= /- C'^(jl,---,in)V'r---C" (31) 

JMa.n 



with the condition with relates, for given p, the indices to the genus g of the 
surface 

n 

(p+l)(25-2 + n) =^(pm,+j, + l). (32) 

1=1 

The cohomology class C^(ji, • • • ,jn) is Witten's class (top Chern class), and 
ipi is the first Chern class. Witten conjectured that F is the string solution of 
the p-th Gelfand-Dikii hierarchy [2, 19]. The partition function Z of (27) is a 
model for the (p,l)-quantum gravity, which is equivalent to the (p-l)-matrix 
model. 

By retaining only the leading term of the limit n — > limit, we restrict 
ourselves to surfaces with one marked point. Then the expansion of (29) is 
simply given by 

lim„_^0--^ = J2 < '^rn,j >g tm,j- (33) 

In the leading order of the limit n — > 0, the matrix A is replaced by a scalar 
A = A • 7. 

After a simple rescaling of the Gaussian weight, we obtain, for instance 
in the case of p=3, 

Z^yJ ^^^^P[^t^B^ + ^trAB^ - itrA^B^] (34) 

Zo^ j d5exp[-^trA2B2] (35) 
In this p=3 case, we obtain from (21), 

1 r(2±i) , , 

< (TW W) ^''^ 

where j = for (7 = 1, 4, 7, 10, • ■ ■ and j = 1 for g — 3, 6, 9, • • •. For g — 
2, 5, 8, ■ ■ ■, the intersection numbers are zero. 

When the genus g is equal to three, the above formula gives < rg,! >g=3= 
^-^^^33,3 which agrees with the value obtained earlier by Shadrin [20]. In the 
calculation of the intersection numbers for p=3, in addition to < (trB^)P >and 
< (trB^)*! > given hereabove in (22) and (26) , one needs to compute mixed 
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averages of the type < (trB"^)'^(trB^)' >. Such averages are indeed required 
if one deals witli tlie expansion of (34), but they are also contained in the 
explicit formula (21) for the generating function. 

For general p, from (21), we have for g=l, < Ti o >g=i= To derive 
this result, we simply need to compute < trB^ >g=i= 1 and < (trB^)^ >g=i= 
3 with ti o = —n/{p\^^^). In the case p=4, we have < rg.o >g=3= g^- 

Thus we find that the expression of (21) for the n — > limit of t/(si, ■ ■ ■ , s^) 
is a generating function for the intersection numbers for the moduli space of 
one-marked point, p spin, curves. This provides thus a simple algorithm for 
obtaining these numbers. 

We have restricted this article to surfaces with one marked point. To go 
beyond one marked point, the analysis of higher orders in n is required. For 
two marked points, we need the order n^. However for genus one, we have 
two different kinds of intersection numbers, < tqT2 > and < rf > (here p=2). 
These two terms are distinguished by coupling to different combinations of 
the ti parameters, namely to^2 and t\ ; those combinations have different A 
dependence (30). Thus we need a matrix A which is no longer a multiple of 
the identity. If we go to higher orders in n with simply for A a multiple of 
the identity matrix, we obtain the values of the sum < T0T2 > + < rf >= j^, 
instead of that of the individual terms. We hope to be able to discuss higher 
marked point by extending the present approach to such cases. 

5 Conclusion 

The present article, using both duality and replica, provides exact formulae 
for the intersection numbers of the moduli of curves on Riemann surfaces 
with one marked point and p-spin. The method may be extended, at least for 
low genera, to a higher number of marked points. It provides a relationship 
between Kontsevich's Airy matrix model and Okounkov's work of the edge 
of Wigner's semi-circle. 
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Appendix A : Gaussian averages in the presence of an external 
matrix source 

For the sake of completeness we reproduce here the main steps of the 
derivation given in [7]. The probability distribution for Hermitian n x n 
matrices is thus 

Zj 

Let us first compute the one point function 

IJa{s) = ^(tre^^)- (38) 

Using the Itzykson-Zuber formula [11, 12] to integrate out the unitary degrees 
of freedom one obtains : 

-1 -| „ n n 

^^^'^ ^ n-zk^) I ^ dhHB)e^^^^-^''^^''^^ E e^"' (39) 

in which Z is a normalization constant fixed such as Ua{0) — 1 and A(S) is 
the VanderMonde determinant of the eigenvalues A(S) = Yli<:j{bi — bj). We 
now use repeatedly the trivial identity 

/ n (i6iA(B)eSi (-^'''+'^^''') = C„eSi A(^) (40) 
1 

in which Cn is a simple constant, with replaced by dj = + s6i^a in which 
a takes the successive values 1 to n. This gives 



Ua{s) = -e^ E 



^ lei V e^««/^ TT """""^^^ (41) 
which may be replaced by the contour integral around the n points (oi, • • • , a„) 

Note that this is an exact formula ; simplifies of course for the pure 
Gaussian case aa = 0, although the derivation does require a source matrix 
in the intemediate steps. 

The derivation for the k-point functions follows exactly the same lines : 
the identity (40) requires now an d which carries k-indices. The result of the 
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calculation is similar : 

C/A(si,---,Sfc) = ^(tre^i^---tre^>^^) 



ns 



{Uj - Uj + Sj - Sj){Ui - Uj) 
i<j {Ui - Uj + Si){Ui - Uj - Sj) 
k n 

xnn(i+^)- (43) 

i=l a=l "a 

One can simplify the expression by noticing the Cauchy determinant identity 

det = (_l)n(n-l)/2 n^<J{x^ - Xj){y. - y,) ^^^^ 

Xi Vj Y[i,j{Xi Vj) 

with Xi = Ui + Si,yi = Ui, namely 

1 ^ 1 TT {Ui - Uj + Si - Sj){Ui - Uj) 

Ui + Si- Uj Si • • • Sfe {Ui - Uj + Si){Ui - Uj - Sj) ' 

In the n = limit we need to consider only the the connected part of Ua (the 
disconnected ones vanishing with n, as is obvious from Feynman diagrams or 
explicit formulae)) which reads, which correspond to connected permutations 
in the expansion of the determinant. Taking a vanishing source (a^ = 0) the 
n limit gives immediately (19). 
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Appendix B : Duality 

Let us consider the Gaussian average, with a matrix source A, of the 
product of characteristic determinants of the N x N Hermitian random ma- 
trices 

1 ^ 

Ffe(Ai, Afe) = ^ < n det(Aa - M) >a,m 

1 k 

= fdMll det{Xi ■ I - M)e-5'^(^-^)' (46) 

in which ^ J dMexp{-l/2TryP) = 2^/2(71)^'/^ ^j^g unitary invariance 
of the measure allows us to assume, without loss of generality, that A is 
a diagonal matrix with eigenvalues (ai, ■ • • , a^). We introduce now N x k 
complex Grassmanian variables {i = 1,2, - ■ ■ ,N and a = 1,2, ■ ■ ■ , k) {i^l, V'f ) 
{i — 1,2, - • ■ ,N and a — 1,2, - • • ,k) with the normalization 

( ) = ( ? ) (*^' 

Then one may write 

k 



n det(A, • I - M) = / n di^td'^t^^''''=' ELi#(^«-^)^..^" (48) 
i=i 

and perform the Gaussian average 

^ / ciMe-^™'+™^ = e^^^l (49) 



We deal here with a matrix X given by 



k 

X,,^a,5,,-Y.r^^l^^ (50) 



and thus 

N k 



TrX^ = TrA^ - 2 5: ^ ap^X - E 7a,/J7 (51) 

p=l a=l a,/3=l, •••,]£ 



where 

AT 



1=1 

To make the notations more transparent let us denote by 'Tr' the trace on 
the initial space ol N x N matrices and by 'tr' the trace in the new space of 
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k X k matrices : the last term in (51) for TrX^ is thus — tr7^, the minus sign 
being due to the anticommuting nature of the psi's. We then introduce an 
auxihary Hermitian k x k matrix /3 so that 

Zk J 

This leads us to the representation 

x/d,e-4-^--<-.^.e-ELEL...#.f (54) 

in which A is the k x k diagonal matrix with eigenvalues (Ai, ■ ■ ■ , A^). One 
can now easily perform the integration over the Grassmanian variables since 
the exponent is quadratic in those variables. This yields the announced dual 
representation 

Ffe(Ai, A,) = ^ /' dPe--^''"" n det[(A, - a,)5,,, + (55) 

in which the matrices are now k x k. 

Remark : this duality takes a more symmetric form if we chose a pure 
imaginary diagonal matrix A, and if we shift the matrix (3 by iA in the final 
formula (55). Then the duality reads 



1 k 

fdMll det{Xi ■ I - M)e-^^(^-^^)' = 
{-if- / dp e-^*^(^+^^)^ n det(a,V - Pp,<r) (56) 



This duality may be extended to a two-matrix model when the measure 
is an exponential of a quadratic form in the two matrices M and M', namely 
exp {aTrM^ + bTrM'^ + cTrMM'). 
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Appendix C: Contour integrals in the n = limit 

We return to the contour integral (19) which gives the n = limit of 
the k-point function. Let us consider for simplicity the case of the two-point 
function. We are then dealing with the sum of two integrals : 



.2 , .2 



lim Uisi, S2) = -e^ / n ^ei:i("'^^A)[log (1 + ^) + log (1 + ^)] 
n^o J 2^7^ Ul U2 

xdet . (57) 

The continuation to n = requires to take contours in the Ui and U2 planes 
which circle around the respective cuts [— si,0] and [— S2,0]. We have to 
choose some well-defined contours on the two variables before we can write 
the integral (57) as a sum of two integrals. For instance we choose and 
integral over a large contour in the U2 plane , and close to the cut in the Ui 
plane. The disconnected term l/siS2 of the determinant gives a vanishing 
contribution. Thus 



.2 , „2 



U{si, S2) = / n |^eE?K-A)[iog (1 + ^) + log (1 + ^)] 

J {- ZlTT Ul U2 

X- i (58) 

[Ul + Si- U2){U2 + S2 - Ul) 

The second integral, the one which involves log (1-1-;^), vanishes for the 
choice of contours that we have made, since we can integrate over Ui first 
and there are no poles inside the countour. For the first part, that with 
log (1 -|- ^), we integrate over U2 first, pick up the two poles and find 

U{si,S2) ^ e-^-ie-'"/^ - e'^'^/^) / ^e^^^^/Mog (1 + ^) (59) 

a J ZZTT Ul 

with a — si + S2- Since 

/ ^gnia/A (1 + fl) ^ _ ^^gWA ^ _ g"^!'^/^) (60) 

J 2m Ul Jsi a 

we end up with 

t^(5i,S2) = 4— sinh^sinh^ (61) 

This calculation may be repeated for the k-point function, although the 
combinatorics becomes heavy. This yields the final result 

lim?7(si,S2,---,Sfe) = — ]j2sinh^ (62) 
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where a = si + S2 + ■ ■ ■ + Sk- 

It is straightforward to generahze this formula to the case of a non-zero 
external matrix source A with eigenvalues (oi, • • • , a„). We assume that the 
density of eigenvalues 

1 " 

p{x) ^ -y25{x - Qi) (63) 

^1 

has a finite limit po{x)'when n goes to zero. Up to now we have dealt with 
Po{x) = 5{x), but we could take other examples such as n/2 eigenvalues equal 
to +a and n/2 equal to —a. Then we would deal with po{x) = ^[S{x — a) + 
5{x + a)]. The derivation goes through now in the same way : we write 

which may be expanded when n goes to zero provided, as we have assumed, 
that p{x) has a limit. The result is then simply 

limt/A(si,S2,---,Sfe) = lini;7(si,S2,---,Sfc) \ dxpo{x)e'"'^^ (65) 

n— »0 n— >0 J 

in which the first factor U is given again by (62). 
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